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O , 

We generalize the concept of the geometric phase in stochastic kinetics to a noncyclic 
evolution. Its application is demonstrated on kinetics of the Michaelis-Menten reaction. It is 
shown that the nonperiodic geometric phase is responsible for the correction to the Michaelis- 
£<\ i Menten law when parameters, such as a substrate concentration, are changing with time. 

We apply these ideas to a model of chemical reactions in a bacterial culture of a growing size, 
where the geometric correction qualitatively changes the outcome of the reaction kinetics. 



I. INTRODUCTION 



Biochemical reactions are typically characterized in stationary in vitro environments with the 



hope that their measured properties will hold in vivo. There are clearly many important physio- 
logical reasons why this extrapolation may fail. In this article, we focus on one particular reason 
that has little to do with the physiology, but rather derives from the fact that rates of complex 
chemical reactions may have non-trivial corrections due to slow, adiabatic drift of (internal) kinetic 
parameters of the system^. 

The class of phenomena we study is related to the celebrated Berry's phase in driven quantum 

•rH . 

^ ■ mechanical systems^, which predicted a contribution to the phase of an adiabatically changing 

wave function in the form of a integral over the parameter trajectory. Since the original Berry's 
discovery a number of its generalizations were proposed, e.g., to nonabelian and nonadiabatic 
regimes. Similar geometric phases were also found in other fields, for example, in dissipative 
dynamics 3 - - — . 

Recently, new geometric phases where studied in the domain of purely classical stochastic 
kinetics^-i^. They were shown to be responsible for the stochastic pump and other ratchet-like 
effects, and thus they are of clear importance for the theory of chemical enzymes, and specifically 
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molecular motors operating in strongly stochastic environmen t 13 1 14 . This finding raises possibilities 
of various generalizations of the geometric phase. For example, recently its nonadiabatic counter- 
part was introduced in Ref.— , and it was shown to be responsible for a non-adiabatic current 
contribution that has no analog under stationary conditions. 

In this Letter, we study another generalization of the geometric phase in stochastic kinetics, 
namely to a nonperiodic evolution in the parameter space. While its quantum and optical versions 
were explored in a series of studies^^i, their role is still largely unclear. In this work, we show that 
the gauge invariant noncyclic geometric phase in stochastic kinetics can be unambiguously defined, 
and that it can be naturally interpreted as being responsible for the leading nonadiabatic correction 
in the expression for stochastic fluxes, which can qualitatively change kinetics of a chemical reaction. 

II. GENERATING FUNCTION FOR THE MICHAELIS-MENTEN REACTION 

The Michaelis-Menten (MM) reaction^ is the most fundamental and the simplest enzymatic bio- 
chemical process. It describes a catalytic conversion of one type of molecules, called the substrate, 
into another type, called the product, via an intermediate reaction with an enzyme. Schematically 
the MM reaction can be represented as 

fcjna fc 2 

S + E^SE^E + P, (1) 

k 2_ 2 n P 

where S and P denote substrate and product respectively, s and p stand for their concentrations, 
and E is the enzyme molecule. S and P interact via creating a complex SE which is unstable 
and dissociates either back into E and S or forward into E and P. In the simplest version of the 
MM mechanism, enzymes catalyze the process but are not modified in any reactions. However, 
generalizations are certainly possible 26 . 

In their 1913 article 25 , Michaelis and Menten considered a strongly nonequilibrium situation, 
neglecting the backwards E + P association, which can be done for n p <C k\n s /k-2- However, here 
we keep this reaction for generality. If the number of S and P molecules is much larger than that 
of the enzymes, the latter have to perform many substrate conversions each in order to change 
S and P concentrations noticeably. This is traditionally used to simplify the reaction kinetics 
since one can assert that enzymes operate in a quasi steady state at current substrate and product 
concentrations. 

Stochastic kinetics of the conversion of S into P is conveniently described by the moments 
generating function Z{x,t) (mgf) and the cumulants generating function S(x,t) (cgf) defined 
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as 7.27.28 



Z( X ,t) = e s ^= J2 P ne mX > ( 2 ) 

n=— oo 

where P n is the probability to find net n product molecules generated during the observation time 
t (back conversion is counted with the negative sign). For a small number of enzymes, they can 
be considered statistically independent over short periods of time, and the cgfs are additive. Thus 
we will restrict our study only to the case of a single enzyme without loss of generality. 

It is convenient to introduce additional generating functions Ue = Yl^L-oo PnE^ nx and Use = 
Y^=-oo PnSE^ mx , where P n E and P n SE are the probabilities that, at a given time, the net number 
of generated product molecules is n, and the enzyme is in the unbound/bound state. Then the 
master equation for the entire process is 

jgPnE = -(hn s + k^ 2 n p )P nE + k-iP n SE + klP(n-l)SE, 
di P nSE = -0-1 + k 2 )P n SE + k\U & P nE + k_ 2 n p P( n+1)E - 

Multiplying ([3]) by e txn and summing over n we find the equation for the generating functions: 



H(x,t)\ , (4) 





where 

/ kin s + k^ 2 n p -k-i - k 2 e %x 

H(x,t)=\ |. (5) 

V —k\n s - k- 2 n p e %x A;_i + k 2 

If we set n = at initial moment t = 0, then the initial conditions for are Ue{P = 0) = pe(0), 
and UsE^t = 0) = Pse(0), where pe(0) and pse(0) are probabilities that the enzyme is free/bound, 
respectively. Additionally, note that Z(x,t) = Ue(x^) + Use(x^)- Thus the formal solution for 
the mgf ([2]) can be expressed as an average of the evolution operator 

Z(x,t) = (l|f(e-J?^*)*)|p(0)>, (6) 

where (1| = (1, 1), \p(0)) = {pe(0),Pse{0)) T , and T is the time-ordering operator. 

Before we proceed with the case where parameters are time dependent, it is instructive to look 
first at the stationary regime. To simplify ©, one can find normalized left and right eigenvectors 
(uq/i\, \uq/i) and corresponding eigenvalues e /i of the operator H{x), where indices and 1 
correspond to the two eigenvalues with the smallest and the largest real parts, respectively. There 
is one left and one right eigenvectors for each eigenvalue. 
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Every vector, such as \p(0)) can be expressed as a sum of eigenvectors of H(x), for example, 

|p(0)> = («o|p(0))|«o> + («i|p(0))|«i>, (7) 

where we define (a|/3) = ol\P\ + a 2 /3 2 to be a standard scalar product of two vectors. Substituting 
d?]) into ©, for the time-independent Hamiltonian we find the steady state mgf, 

Z st (x,t) = e - e o(x)*+ln((l|«o){«o|p(0))) + e -ei(x)t+ln((l|ui>{uib(0))) ) ^ 

At time scales t S> max[l/A;_i, 1/&2; l/{k\n s ), l/(fc_2n p )], the second term in ([8]) is exponentially 
suppressed in comparison to the first, and the expression for the mgf simplifies to 

Z st (x,t) « e — e o(x)*+in({i|«o)<«ob(o)))_ ^ 

Terms analogous to — eo(x)t in © have been studied previously"^. The second term is less 
threaded: this is the boundary term that does not grow with time and depends on the initial 
conditions and the averaging over the final states of the enzyme. One can disregard it in comparison 
to the first contribution when t — > 00. However, we note that its relative effect decays as 1/t, i.e., 
not exponentially. We will keep the boundary term in the following discussion because it will play 
an important role to restore the gauge invariance of the nonperiodic geometric phase. 

At the first look, the boundary term leads to a contradictory result after setting t — > 0, i.e at 
the initial moment of the evolution. In this limit, the boundary term does not disappear, namely 

SbndU = In«l|«o(0)><«o(0)b(0)» + 0. (10) 

However, we expect Sbndlt=o to be zero, since n\ t=0 = 0, so the mgf should be identically equal to 
unity. The apparent contradiction is resolved by noting that @ was derived assuming t — > 00, and 
it is simply an invalid approximation for t = 0. In other words, the boundary term is responsible 
for the initial fast relaxation to the stationary regime. For more insight, one can calculate the 
contribution of the boundary term to the average number of generated product molecules. Using 
the normalization condition Pse(0) = 1 — Pe(0) one can find 

9 <Sbndlt = o 



"-bnd = ~« 



_ (fc 2 + fc_ 2 n p )(fc 2 + fc_i - Kp E (0)) 

K 2 ' ' ' 



where K = k—\ + A; 2 + k\n s + /c_ 2 n p . If one assumes that the initial probability pe(Q) for the 
enzyme to be free is at the equilibrium value pe(0) = (^2 + k-\)/K, then (fTTj) produces nbnd = 0, 
as expected. To confirm this, one can also derive by a standard master equation approach. 
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That is, calculating the average number of new product molecules Wbnd(^); one would find that, 
after a sufficiently long time, 

n(t) = n bnd + - p t. (12) 

The second term in Eq. (|12|) is the average number of the product molecules produced during time 
t at a steady state. It is the standard prediction of the reversible MM theory, and the first term is 
a correction, which is nonzero when the initial state of enzymes is not the same as its steady state. 



III. NONCYCLIC GEOMETRIC PHASE IN STOCHASTIC KINETICS 

Assume now that there are several slowly time-dependent parameters in the model. We will 
group them in a vector A. In the case of the MM process, one can view these time-dependent pa- 
rameters as concentrations of the substrate and the product, A = (n s ,n p ). However, the discussion 
in this section is completely general. 

Following Ref*£ we partition the time into small intervals, over which kinetic rates can be 
considered almost constant. We insert the resolution of the identity operator, 1 = \uo(t))(uo(t)\ + 
\ui(t))(u\(t)\, in ([f)D after every such an interval. One can find then that the boundary term 
becomes Sbnd = ^ n ((M u o{t)){uo{0)\p(0))). Importantly, it is no longer gauge invariant, i.e., it is 
sensitive to the redefinition of eigenstates of the Hamiltonian ([5]) such as \uq) —> e^'|uo) and 
(uo\ —> (uo|e -< ^ A ). Therefore, taken alone, it has no direct physical meaning. 

It will be convenient to rewrite the boundary term as a sum of a gauge invariant part and a 
term that is an integral from a pure derivative, i.e., 

Sbnd = S hnd \ t=0 + J P • dX, P = d x ln(l|« ), (13) 

where c is the contour in the space of the variable parameters. By analogy with Refi^, and 
including the boundary contribution (|13p . the mgf in the quasi steady state limit can be written 
as an exponent of a sum of two terms, 

Z(x) = e Sgcom( -^ +Sqat( - x \ (14) 

where 

Sqst(x) = - € (x,t')dt' + Sbnd| t =o ( 15 ) 

Jo 

is the quasistationary part of the generating function averaged over time. This is the part that 
morphs into the steady state result © for fixed values of all parameters. 
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The other term in ([T 



S g eom = /[P(A) - A (A)] • dX, A (A) = (u \d x u ), (16) 

is the geometric phase contribution responsible for additional reaction events. A is called the Berry 
connection. S ge0 m has no analog in the strict steady state regime. 

Note that, unlike in Ref^, we do not assume a periodic evolution of parameters. Therefore, 
the term involving the integral of the Berry connection A over the path in the parameter space, 
— J c A(A) • dX, is generally not gauge invariant. However, one can easily check that the non-gauge- 
invariant contribution due to the boundary term exactly cancels the non-gauge-invariant part of 
the contour integral from A. 

We further mention that the definition (|16p differs somewhat from those used for the non-cyclic 
geometric phase in quantum mechanics. For example, Refs i 16 ' 17 define the noncyclic geometric 
phase as 7 gp = J* c [A(A) — P(A)] • dX, where P = — Im ( ^1)a°(o)^m?a^ ) ■ ^ n * ne P resen t context, the 



meaning of such definition is unclear, while the geometric phase defined in (|16p is derived directly 
from the exact representation of the mgf. 

Since P is a pure gauge, it is important only when looking at an evolution along an open path 
in the parameter space. If the parameter vector A returns to its initial value at the end of the 
evolution, the expression (fl6|) becomes equivalent to the full-period geometric phase defined in 
RefJ. 

IV. CORRECTIONS TO MICHAELIS-MENTEN LAW 

Consider now the average product creation rate in the MM system under the slow parameter 
evolution. The average number of new product molecules is (n(t)) = —i(dZ(x,t)/dx) x= Q- There- 
fore, just like the full cgf, the average rate of the product production (J) = d(n(t))/dt can be 
written as a sum of the quasistationary J qs t and the geometric J ge om contributions 

de (x,t) 



I T\ — T _i_ 7 _ A dSgeom 
\J i — •-'geom ~r ■-'qst 



dx 



(17) 



dt dx 

The geometric phase is time-dependent only via the time-dependence of the parameter vector A. 
In the case of MM reaction with time-dependent concentrations re s and n p , the time derivative of 
the first term in (fT?]) can be expressed as d/dt — > (dn s /dt)d/dn s + (dn s /dt)d/dn p . Substituting 
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the eigenvectors and eigenvalues of H(x, A) into (fT7|) . we find 

_ (km s (t))k 2 - (fc_ 2 n p (t))A;_i 
J qst - ^ , 

</geom = -(h + k-l) ^3 • (19) 

One can recognize J qs t as the average current for a steady state with fixed values of parameters. In 
fact, (|18p is what is known as the Michaelis-Menten law. However, our results show that this law 
is not exact when concentrations of the substrate and the product have their own time-dependent 
evolution. The geometric contribution is the first correction to the Michaelis-Menten kinetics that 
becomes nonzero when the substrate/product concentrations change with time. Specifically, in the 
most frequent case n p ~ 0, the average rate of the coarse grained MM reaction per one enzyme 
becomes 

, T \ _ k 2 n s k 2 k 1 h s (t) 

' J >"^TfeTET- '^-' W + fe + fc-!) 3 - (20) 

That is, even in this case, the time-dependence of the substrate concentration introduces corrections 
to the reaction rate. 

It is possible to understand the result (|19p with a simpler approach, which, unfortunately, is 
hard to generalize for higher current cumulants to demonstrate the geometric nature of the effect 
for all of them. The probability pe of the enzyme to be unbound evolves according to the master 
equation 

j t PE = -[hn s (t) + k- 2 n p (t)}p E + {k 2 + fc_i)(l - p E ), (21) 



with the solution 

Pe 



(t) = (k 2 + k-l) f\~ It\^Mr)+k- 2 n p (r)+k 2+ k^]dr dt ^ ^ 

Jo 



The lower limit in this integral is not important because we work in the adiabatic approximation, 
which means that the information about the initial state is quickly forgotten. The exponent of the 
integral over r in (|22p is then 



- / t * [k 1 n B (T)+k-2n p (r)+k2+k-i]dT 

e -[fein s (t)+fc_ 2 n p (t)+fc 2 +A.-i](t-ti) A + k in s {t) + k^ 2 n p (t) ^ _ ^ 2 j ^ 



Performing the remaining integration we find the expression for the probability of the enzyme to 
be unbound: 

^2 + ^-1 , (h + k-i)(kih s (t) + k- 2 hJt)) 
PE*^^ + ^3 • ( 24 ) 
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From (|24|) . one can calculate the average reaction rate and check that indeed, it is the sum of the 
quasi-stationary and the geometric components determined in (|18|) and (|19p . 



J(t) = (1 -pE(t))k 2 - PE{t)k- 2 n v (t) = J qst + Jgeom- (25) 



V. GEOMETRIC RATE CORRECTIONS IN A GROWING CELL CULTURE 



The geometric correction (|19p is generally much smaller than the main contribution (|18p if the 
number of the enzymes is much smaller than that of the substrates and the products. However, 
this small correction has very different properties, and can change a system behavior qualitatively 
under special conditions. 

The quasi-steady state contribution to the kinetic rate in ()18p can be vanishing due to a sym- 
metry relation, such as the detailed balance condition, which guaranties that all chemical fluxes 
at the thermodynamic equilibrium state are zero on average. Thus, if a system is slowly driven 
externally so that it always remains close to the thermodynamic equilibrium, the quasi-steady state 
approximation will predict zero average product creation. In contrast, the geometric contribution 
does not have to remain zero, and it will result in a qualitatively novel effect. 

To show this, consider the MM reaction with concentrations of substrate and product n s and n p 
large and treated deterministically. Let us suppose that the system is initially in an equilibrium, 

hk 2 n s (0) = fc_i/c_ 2 n p (0). (26) 

Now suppose that this process happens inside a living cell that grows and divides in its usual cycle. 
Then the substrate/product molecules in a single cell are diluted by N(t), the number of cells in 
the descendant colony: 

n s (t)=n s (0)^||, n p (t)=n p (0)^||, (27) 

Since the ratio n s (t)/n p (t) is not affected by this time dependent dilution, the system remains 
near equilibrium, and the quasi-steady state reaction rate remains zero. Thus the average number 
of new product molecules, produced by a single enzyme is completely determined by the geometric 
part of the rate ([!§]). 

"OO 

n = I dv 



»(0) 



ih + k - l) KHv) 

hk 2 n B (0) 



(k 2 +k-i)(k-i +fc 1 n B (0)) 



• (28) 
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On the one hand, this effect is very small: the average number of new product molecules per 
one enzyme is a fraction of unity, which compares to a large number of already existing substrate 
and product molecules. On the other hand, the geometric contribution qualitatively changes the 
result, predicting on average nonzero amount of new product molecules, which is not expected from 
the standard MM treatment. If the number of the enzymes in the culture is proportional to the 
number of cells, and hence grows with time as N(t), this effect will eventually become observable. 

The result (f25j) would be valid only if we could treat concentrations as parameters, changing 
only due to the external volume growth. In a closed system chemical fluxes eventually should be 
compensated by the reverse fluxes due to the violation of the steady state condition (f26|) . Thus 
the geometric flux should be possible to detect by measuring the deviation of the ratio n s /n p from 
the equilibrium value. 

Considering intermediate stages of the culture growth, one can notice that the number of newly 
produced molecules depends only on the initial and the final cell numbers: that is, the average 
number of produced proteins depends on the current state of the system, but not on how it got 
there or where it's going from there. This can be utilized by living organisms in order to control 
some processes depending on the stage of cell's life cycle. Although this effect is very small, it 
should be interesting to explore its detectability in vivo and employ it in artificial biochemical 
circuits design. 

VI. DISCUSSION 

In this article, we generalized the notion of the geometric phase in evolution of the mgf to 
nonperiodic time-dependent processes. For this, the contour integral of the Berry connection had 
to be supplemented by an extra term restoring the gauge invariance of the geometric contribution 
to the cgf. This term originates from the boundary contribution responsible for proper description 
of the initial and final moments of the measurement. For nonequilibrium initial conditions, the 
boundary terms are responsible for the initial fast relaxation to the enzymatic quasi-steady state. 
That is, although our approach is adiabatic, it also rigorously captures initial fast relaxation effects. 

Our non-periodic geometric phase is different from the ones often encountered in quantum 
mechanical applications. Its uniqueness follows from the existence of a special gauge that should 
be imposed in order to describe stochastic kinetics correctly. 

We showed that the phase is responsible for nonadiabatic corrections to the standard Michaelis- 
Menten approximation. Such corrections are usually small in comparison to the quasi-steady 
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state predictions. However, they explicitly break time-reversal symmetries and, therefore, can 
produce a qualitatively different result when a chemical system is driven closely to a thermodynamic 
equilibrium, as in the cell culture growth model that we discussed. 

It is unclear whether this effect is of biological relevance. However, we note that we studied 
only the simplest of its realizations. The introduced non-periodic geometric phase is completely 
general and should appear practically in any interacting chemical system driven by external fields. 
Other interesting examples will surely emerge with time. We expect the greatest opportunities for 
biological relevance in the domain of molecular motors, where geometric effects play an important 
role as is^. 

It would also be interesting to understand if the nonperiodic geometric phase is related to the 
existence of fluctuation theorems^. Indeed, instead of chemical fluxes, it is possible to use the 
same formalism to count work or dissipated energy in a driven stochastic system. The absence of 
anholonomies, such as cyclic geometric phases may indicate the existence of fluctuation relations 
since then the counting statistics depends only on initial and final values of external parameters, at 
least in the adiabatic limit. Generalizations of our approach to a nonadiabatic evolution should also 
be possible since similar generalizations simultaneously to a noncyclic and nonadiabatic evolution 
in quantum mechanics exist^S. 
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